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ABSTRACT 


Calculations of the steady-state photoelectron energy and angular 
distribution in the altitude region between 120 and 1000 km are presented. 
The distribution is found to be isotropic at all altitudes below 250 km, 
while above this altitude anisotropies in both pitch angle and energy are 
found. The isotropy found in the angular distribution belov/ 250 km Implies 
that photoelectron transport below 250 km is insignificant, while the 
angular anisotropy found above this altitude Implies a net photoelectron 
current in the upward direction. The energy anisotropy above 500 km arises 
from the selective backscattering of the low energy photoelectron popula- 
tion of the upward flux component by Coulomb collisions v/ith the ambient 
ions. The total photoelectron flux attains its maximum value befeveen about 
40 and 70 km above the altitude at which the photoelectron production rate 
is maximum. The displacement of the maximum of the equilibrium flux is 
attributed to an increasing (with altitude) photoelectron lifetime. Photo- 
electrons at altitudes above that where the flux is maximum are on the 

average more energetic than those belov/ that altitude. The flux of photo- 

8 

electrons escaping to the protonosphere at dawn was found to be 2,6 x 10 

- 2-1 8-2 
cm'’^ sec” , while the escaping flux at noon was found to be 1.5 x 10 cm 

-1 9-2 

sec~ . The corresponding escaping energy fluxes are: 4.4 x 10 eV cm 

sec’^ and 2.7 x 10^ eV cm“^ sec"\ 
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1. INTRODUCTION 

To study dayglow emissions, electron and ion temperatures, predawn 
airglow intensity and electron temperature enhancements, exospheric 
neutral particle temperatures , etc. one needs to know the equilibrium 
energy spectrum of the photoeVectrons (energetic electrons) that consti- 
tute the energy source of these phenomena. Several theoretical studies 
(Hoegy et , 1965; Nisbet, 1968; Henry and McElroy, 1968; Dalgarno 

, 1969; Rees 6t , 1969; Banks and Nagy, 1970; Nagy and Banks, 
1970; Takayanagy and Itikawa, 1970; Cicerone and Bowhill, 1971a; 

Dalgarno and Lejeune, 1971) have progressively improved our knowledge 
of the energy spectrum of photoelectrons in the ionosphere. In parallel,., 
experimental studi.es (Shea et , 1968; Yngvesson and Perkins, 1968; 

Rao and Maier, 1970; Heikkila, 1970; Evans and Gastman, 1970; Doering 
^ , 1970; Cicerone and Bowhill, 1971b; Knudsen and Sharp, 1972; 

Galperin ^ al_. a 1972) have provided valuable data for comparison with 
the theoretical results. 

In a recent review Cicerone et al . (1973) have found that at high 
altitudes differences of about a factor of two are found among photo- 
electron fluxes calculated by different theoretical methods. Although, 
these differences can be eliminated within the existing theoretical 
framework, experimental techniques (Doering ^£!_. , 1973; Hays and 
Sharp, 1973) have improved substantially in resolution and accuracy 
to warrant more detailed theoretical calculations than hitherto under- 
taken. In a previous paper (Mantas 1974) a detailed theoretical 
formulation of the problem of photoelectron thermal Izati on and transport 
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in the Ionosphere was presented. In the present paper we use that 
'theory to calculate the steady-state pitch angle and energy distribution 
of photoelectrons for two ionospheric models appropriate for conditions 
prevailing over Arecibo (18.4°N, 67° W) on December 24, 1968 at dawn 
(0703 hours, solar zenith angle x ~ 90°), and on June 26, 1968 at noon 
(1225 hours, x : 0°). 

The boundary conditions for the photoelectron thermal izati on and 
transport equation have been discussed by Mantas (1974). In the 
present calculations we have assumed (upper boundary condition) that there 
is no photoelectron flux incident upon the ionosphere at 1000 km. At 
120 km (lower boundary) a partially reflecting boundary, with a reflec- 
tion coefficient of 0.6, was assumed. It was found, however, that the 
lower boundary condition has no significant effect on the solution because 
at this altitude transport is insignificant. The contribution of a 
non- zero (downward flux at the upper boundary can be calculated separately 
and added to the present calculations. 
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2. IONOSPHERIC MODELS, PHOTOELECTRON PRODUCTION RATES. AND CROSS SECTIONS 

2.1 Neutral Atmospheres. 

The neutral particle densities used in these calculations have been 
obtained from the Bates (1959) model atmosphere as elaborated by Stein 
and Walker (1965) and Walker (1965). The neutral gas densities at 120 km 
and the other parameters of the Bates model have been obtained from a 
combination of theoretical considerations and temperature data from 
Arecibo at the indicated times (R. J. Cicerone, Private communication, 1972) 
The neutral atmosphere exospheric temperature for the sunrise model is 
884 Kelvin, while that of the noon model is 1114 Kelvin. 

2.2 Charged Particle Densities. 

The electron density profiles have been obtained from Thomson scatter 
observations over Arecibo at the indicated times (R. J. cicerone. Private 
communication, 1972). The electron and neutral atmosphere densities can 
be found in Mantas (1973) . Since all ions are assumed to be singly 
charged and their mass is much larger than the electron mass, all electron.. 
Ion collisions are characterized by the same collision cross section and 
no distinction is made between different kinds of ions, 

2.3 Primary Photoelectron Energy Spectra. 

Energy spectra of the primary photoelectron production rates, at 
selected altitudes, for the sunrise model are shown in Figure 1. These 
spectra have been calculated by R. S. Stolarski (private communication, 1972) 
using the EUV flux measurements of Hinteregger (1970) and the photo- 
ionization cross sections reviewed by Stolarski and Johnson (1972). 

The pitch angle distribution of the primary photoelectron production rate 



is assumed to be that given by Mariani (1964). 
.2.4 Cross Sections 
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We have used the empirical formula 
q„A. W.5j . 1, 

o — 0 - c ) (c 


proposed by Green and Dutta (1967), to represent the excitation cross 

sections. In equation (T) E is the electron energy, is the excitation 

-14 2 2 

potential of the state j, = 6.51 x 10 cm eV , and A^, gy are 
dimensionless parameters obtained by fitting. (1) to experimentally 
determined cross sections. For the excitation of atomic oxygen we have 
used values for the parameters of (1) suggested by R. s. Stoiarski (private 
conrouni cation, 1972). For molecular nitrogen, molecular oxygen, and 
helium we have used the values given by Stoiarski et ^. . (1967), Watson 
£t (1967), and Jusick (1967), updated to more recent cross 

section measurements by Jobe ^ (1967), Burns ^ al_. (1969), McConkey 

and Simpson (1969), Simpson and McConkey (1969), Stanton and St. John 
(1969), Brinkman and Trajmar (1970), Sheridan ^ aj_. (1971), McConkey 
and Woolsey (1969), and Skubenich (1968). The probability densities - 
p|(£',E) (see Mantas 1974^ for the energy distribution of the post- 
collision electrons in electron-neutral ionizing collisions have been 
constructed from the differential cross sections given by Khare 
(1969) and Dalgarno and Lejeune (1971). The total ionization cross 
sections have been taken from the review paper of Kieffer and Dunn (1966). 
A listing of the values of the parameters of (1) is given by Mantas (1973). 
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3, RESULTS AND DISCUSSION 

3.1 Photoelectron Flux Pitch Angle and Energy Spectra at Dawni 

In Figures 2a through 2d, the logarithm of the photoelectron flux 
(log^Q 2 tt $(h,E,|i)) is plotted at selected altitudes. For negative 
values of the direction variable (y < 0) the. electrons are moving in 
the downward direction along the magnetic field lines, while for posi- 
tive values they move in the upward direction. Note that the flux shown 
in these figures has been integrated over the azimuthal angle the 
quantity 2n 4>(h,E,y), therefore, represents the flux over a pitch angle 
cosine interval Ay and not over a unit solid angle. 

Intercomparison of these figures reveals variations of the flux with 
each of the independent variables. For altitudes bdlow about 250 km 
4-(h,E,y) is isotropic (independent of y).- For altitudes above 250 km 

more electrons are moving in the upward than in the downward direction 

the anisotropy becoming increasingly larger with increasing altitude. 

It is important to note that most of the flux variation occurs in the 
small interval -0.25 £ y ;< 0.25 or, in terms of the pitch angle^in the 
region between about 75 and 105 degrees. Outside this region the 
dependence of $ upon y is weaker. Above 800 km the angular dependence- 
is much stronger; however, in this region most of the variation is intro- 
duced through the somewhat unrealistic upper boundary condition that no 
downward flux is present at the 1000 km level. Although at certain times 
of the year conditions exist where no flux from the conjugate ionosphere 
is present, some of the upward flux will be scattered downward at altitudes 
higher than 1000 km, and thus constitute a nonzero upper boundary condition 
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An estimate of the reflection coefficient at the upper boundary will be 
^iven later. 

Now we turn our attention upon the examination of the photoelectron 

/ ‘ 

flux energy spectrum. Examination of Figures 2a through 2c reveals con- 
siderable structure of $(h,E,y) with respect to the energy variable. 

For altitudes below about 250 km the local energy variation is much 
smoother than that above this altitude. A comparison of Figure 2a 
with the primary photoelectron spectrum at low altitudes (see Figure 1) 
reveals that the energy spectrum of the steady -state photoelectron flux 
is much smoother than the energy spectrum of the primary production rate. 
The prominent peaks of the primary spectra are absent in the equilibrium 
flux spectra. The wide maximum in the energy region between 20 and 65 eV, 
in the primary spectra, is substantially reduced. The deep and wide 
minimum, in the primary spectra, in the energy about 10 eV,has disappeared 
completely and the flux has a negative slope in the energy interval 
between 10 and 20 eV, while the production spectra have a large positive 
slope in this interval. For altitudes below about 250 km (see Figure 
2a) the steady state photoelectron flux spectra show an increasingly 
larger slope, with decreasing energy below about 40 eV (small local 
variations in the slope are neglected). At about 10 eV the photoelectron 
flux begins to Increase rapidly with decreasing energy. This increase 
is Interrupted at about 5 eV, and begins again, at a much higher rate, 
for energies below about 2,5 eV. 

For altitudes at about 250 km and above, the most prominent featuces 
of the primary photoelectron spectra are also present in the photoelectron 
flux spectra, though less prominent in the latter case. The substantial 
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structure found in the production spectra in the region below about 
'20 eV is absent in the flux spectra* A careful examination of the 
Figures 2b through 2c reveals many interesting variations of the photo- 
electron flux with altitude, pitch angle, and energy. In addition to 
the prominent upward -downward anisotropy that has been mentioned, the 
following important features should be noticed: 

(a) The steep rise in the flux spectrum observed at altitudes 
below 250 km at about 2 eV is absent from the energy spectra 
at about 300 km and above. 

• (b) For altitudes at about 450 km and above, in the energy region 

below about 5 eV^^the upward moving flux decreases with decreasing 
energy the decrease becoming larger with Increasing altitude. 

(c) The features present in the spectra around 25 eV, for a fixed 
direction (v = constant become smoother with Increasing 
altitude above 300 km. 

(d) The smoothness of the spectrum features at a fixed altitude 
increases with decreasing absolute value of the direction 
variable \i. 

These results can be interpreted in physical terms as follows. 

At low altitudes, due to the higher concentration of the neutral atmospheric 
constituents relative to the thermal electron concentration, the photo- 
electron flux energy spectrum is mainly determined by inelastic collisions. 
For altitudes above about 300 km the photoelectron-thermal electron 
interaction becomes increasingly important with increasing altitude. 
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The pitch angle isotropy of the electron flux spectra below about 
250 km shows that electron transport at these altitudes is negligible. 

On the other hand, the increasing pitch angle anisotropy with altitude 
above 300 km shows that electron transport increases in importance as the 
altitude increases . The importance of transport at high altitudes can 
also be deduced from the fact that the electron flux changes little with 
increasing altitude above 400 km, although the photoelectron produc- 
tion rate decreases rapidly. Therefore, to interpret the charac- 
teristics of the electron flux spectra and obtain a qualitative evaluation 
of the validity of the solution, one must consider that most of the 
electrons found at high altitudes are produced in the altitude region 
between 250 and 350 km and are subsequently transported to higher 
altitudes by spiralling around the magnetic field lines. 

The energy lost by an electron in its upward motion through the 
ambient electron gas is a function of the distance it traverses and its 
initial energy. If we' neglect elastic scattering by the ambient ions 
the distance traversed by an electron in reaching certain altitude, 

- say 700 V km, is inversely proportional to the cosine of its pitch angle. 
Although elastic scattering smears the distance traversed by the electrons 
found at the same altitude in the same pitch angle interval, on the 
average one would expect that electrons with large pitch angles, say in 

p 

a small interval around 7Q degrees, would have traversed longer distances. 
In reaching the said altitude, than electrons with smaller pitch angles. 

If the electrons in both pitch angle intervals had originated from the 
same altitude and had the same initial energy, upon reaching a high 
altitude the electrons with large pitch angles would have lost more of 
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their initial energy than those with smaller pitch angles. 

From these considerations one would expect that prominent features 
in the photoelectron production spectra, such as the peaks around 25 eV 
due to the absorption of the He II lines of the solar EUV, if they 
appear in the steady-state electron flux spectra should increase 
in width and shift towards lower energies with increasing altitude and 
pitch angle. The broadening and shift towards lower energies of these 
features as |y| ^ 0 is clearly seen in Figure 2c. 

It should be also observed that at high altitudes the slope 
— $(h,E,y) of the energy spectrum of the upward (y>0) electron flux 
component at the low energy limit (E < SeV) is positive, while that of 
the downward flux (u < 0) is negative. This behavior indicates that the 
upward electron flux component loses low energy electrons at a higher 
rate than it gains, while the opposite holds for the downward flux 
component. Since at high altitudes photoelectrons are losing energy 
mainly to the thermal electron gas and are redistributed in pitch angle 
by Coulomb collisions with the ambient ions, the interplay of these 
two processes must, to a large extent, determine the shape of the spec- 
trum at the low energy limit. The rate at which a fast electron is 
losing energy to the thermal electron gas with distance is, approximately, 
inversely proportional to its energy. And the cross section for Coulomb 
collisions is inversely proportional to the square of the energy. Energy 
loss to the thermal electron gas, therefore, constrains low energy photo- 
electrons from reaching high altitudes and. In addition. Coulomb 
collisions with ions selectively backscatter low energy electrons from 
the upward to the downward flux component. These two processes operate 
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on the downward flux component in the same manner. However, due to the 

larger number of upward moving electrons. Coulomb collisions give rise 

to a net transfer of electrons from the upward to the downward flux 

.component. Since the cross section for Coulomb collisions depends on 

the inverse square of the energy, low energy electrons are preferentially 

backscattered . The net effect' of this 

process is reflected in the shape of the spectra in the low region and, 

as we shall see below, gives rise to an energy anisotropy between the 

upward and downward flux components. 

3.2 Photoelectron Flux Energy Spectra at Dawn 

Integration of 2n $(h,E,u) over p gives the number of electrons per 

unit energy interval crossing a sphere of unit surface in all directions 

per second. This quantity will be referred to as the total photoelectron 

flux energy spectrum and will be written as $^(h,E). Thus 

fl 

$''(h,E) = 27 t $(z,E,w)dji= 4'(h,E) + $'^(h,E) (2) 

j-1 ' 

where n 

r(h,E) = 2 t, ^(h,E,p)dy • (2a) 

.-1 

represents the energy spectrum flux crossing the downward hemisphere, 
and f 1 

**(h,E) = 2it 4(h,E,y)du (2b) 

, 0 

represents the energy spectrum of the flux crossing the upward hemisphere. 
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The total photoelectron flux energy spectrum Is shown in 

Figure 3. It should be noted that this surface is somewhat distorted 
because unequal altitude intervals have been plotted along the depth 
axis as if they were equal. For example <j>^(h,E) has been evaluated at 
9 altitudes in the altitude interval between 200 and 300 km, while in 
the interval between 300 and 400 km it has been evaluated at 7. The 
ratio of the intervals along the depth axis representing the altitude 
intervals from 200 to 300 km and from 300 to 400 km, is 9/7 instead of 
one. This distortion does not affect the accuracy of this figure, it 

o 

should, however, be kept in mind when Figure 3 (and also Figures 5,6 
and 8) is examined. Figures 2a through 2c do not suffer from this dis- 
tortion because, although *{h,E,K) was evaluated at unequal pitch angle 
cosine intervals, interpolation has been used to obtain values of 

$(h,E,u) at equal intervals along the p (depth) axis before plotting^ 

\ 

i 

A careful examination of Figure 3 shows that the sharp peak of the 
spectrum at about 25 eV becomes broader and shifts towards lower energies 
with increasing altitude as discussed previously. The spectrum is also- 
seen to become smoother with increasing altitude above 300 km. 

To obtain a quantitative view of the up-down energy anisotropy we 
evaluate the mean energies; 


E- 
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where: 


$"(h) = 
*‘'(h) = 


t'^(h) = 


j 

i 


100 

^$'(h,E)dE 

100 

4‘^(h,E)dE 

1 

100 

$'^(h,E)dE 


(4a) 

(4b) 

(4c) 


are the total number of electrons crossing a sphere of unit surface, at 
the altitude h per unit time interval, (a) in the downward direction, 

(b) in the upward direction, (c) in both directions. 

The mean energies E^(h) and E^(h) are shown in Figure 4, together 
with the ratio r^(h)/E’(h) and the mean energy of the primary photo- 
electron spectrum (source). From this figure it can be seen that for 
altitudes above about 550 km the upward flux consists of more energetic 
electrons than the downward flux. This is exactly what one would expect 
in accordance with the arguments presented in the previous paragraphs. 

For altitudes between 250 and 550 km the downward flux is composed of 
slightly more energetic (E^ > E^, by about 0.5 eV at most) electrons 
than the upward flux. Below 250 km energy isotropy is seen to 

prevail, as one would expect, due to the large number of collisions at 
low altitudes. The energy anisotropy at high altitudes is probably 
higher than that calculated here, since it is known that at high energies 
the elastic collision cross sections are strongly peaked in the forward 
direction, while we have assumed isotropic scattering for elastic 
collisions independently of the collision energy. The mean energy of 
the total photoelectron flux is found by integrating equation (3c). 

The curve for has not been plotted, since it follows closely the 
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curve for E^. 

« 

It is interesting to compare the mean energy Ej (source) of the 

photoelectron production rate with ^ ~ 'F’. For altitudes above 350 km 

Eg is approximately independent of altitude and has a value of about 18 

eV. Below 350 km rises rapidly with decreasing altitude to a value 

of about 41 eV at 170 km, while decreases with altitude from a value 

of about 17 eV at 1000 km to about 7 eV at 160 km. The large difference 
of about 17 eV at 1000 km to about 7 eV at 160 kra. The large difference 

between and at altitudes below 300 kra indicates that the high energy 

rate as rapid as the rate of production. Above 300 km Inelastic collisions 
cease to play the dominant role, and the equilibrium spectrum is shifted 
only slightly toward the low energy region, since elastic collisions 
with ions absorb no energy and photoelectron-thermal electron energy trans- 
fer is a slow process for photoelectron energies greater than few electron 
volts. 

3.3 Photoelectron Flux Pitch Angle and Energy Spectra at Noon. 

In the present and following subsections we present electron flux 
spectra for noontime conditions. We shall limit the discussion, however, 
only to pointing out what is new or different from the dawn model 
spectra. For details theonterested reader should consult Mantas (1973), 

The main difference in the input par^eter of the two models consists of 
higher neutral and charged particle concentrations in the noon than in the 
dawn model. The F2-layer thickness is larger in the noon than in 
dawn model, with a peak electron density of 10^{cm“^) at about 400 km, 
while the maximum electron density in the dawn model is about 5X10^, 
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and the F^-l ayer peak occurs at about 250 km. Another important dif- 

ference is the higher primary photoelectron production rate, with a 

maximum at a lower altitude, in the noon than the dawn model. The 

3 -3 -1 

maximum primary photoelectron production rate at noon is 3.2X10 (cm sec ) 
at about 150 km, while that at dawn is 1.5 X 10'"(cm sec ) at about 250 km. 

The higher neutral and charged particle densities, and the lower 
altitude of the maximum primary photoelectron production rate in the 
noon model give rise to the following differences between the dawn and 

noon electron flux spectra: 

* 

a) The angular anisotropy at high altitudes at noontime is smaller than 
that at dawn. 

b) . The energy spectrum above 300 km at noon is smoother than the dawn 

spectrum. 

c) The broadening and the shift towards lower energies of the prominent 
features of the spectrum, at fixed altitudes above 400 km, with 
decreasing |u| is larger at noon than at dawn. The same Is also 

true for the broadening and shifting of these features with increasing 
altitude for fixed values of p. 

d) The amplitude of the flux above 400 km at noon is smaller than that 
at dawn, 

A better illustration of the upward-downward electron flux anisotropy 
at high altitudes is obtained by integrating the flux $(h,E,p) over all 
energies. The function 

flOO 

^(h,p) = 2fr $(h,E,p)dE (5) 

‘ 1 

represents the total number of electrons, irrespective of their energy, in 
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the pitch angle cosine internal dp about p. This function is shown in Figure 5. 
It is interesting to note that the pitch angle dependence shown in this figure 
is in qualitative agreement with the single energy analytical results of 

stolarski ;(1972) . The rapid variation of the pitch angle distribution in the 90 

{ 

degree region at high altitudes is clearly seen. It should be also 
observed that outside this region the upward flux component is almost 
isotropic. 

In Figure 6 the total photoelectron flux energy spectrum [$*^(h,E)] 
at noon, defined by (2) is shown. Comparison of Figures 3 and 6 shovis 
that the broad maximum between 200 and 300 km in the high energy region, 
seen in the dawn model, is barely present in the noon model. Another 
interesting characteristic to be noted is the very rapid transition in 
the smoothness of the noon spectrum at about 300 km. In the dawn 
spectrum (Figure 3) no such sharply defined transition altitude is found. 

Further comparison of these figures shows that the broadening and shift 
of the noon spectrum features, with increasing altitude, is larger than 
that for Wwn. A more clear illustration of this effect is given in 

i 

Figure 7 where the energy spectrum of the upward photoelectron flux at 
noon is shown. The prominent peak at 22.5 eV at 304 km is gradually 
broadened and shifted to about 20 eV as the altitude increases to 1000 km. 

The effect of energy absorption by the thermal electron gas and the 
Coulomb scattering by the ions can be seen in the energy dependent attenuation 
of the flux with altitude. For example the flux amplitude in the 
energy region around 70 eV decreases only by about a factor of 2 
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as the altitude increases from 300 to 1000 km, while in the region around 

♦ 

5 eV it decreases by about a factor of 20 over the same altitude interval. 

In Figure 8 the mean energy of the photoelectron flux at noon, as 
function of pitch angle and altitude, defined by 

rlOO 

^0 = Ikiu J/ 

is shown. Below about 220 km the mean energy of the photoelectron flux 
1s independent of altitude and pitch angle and has a value of about 8 eV. 
Between 220 and 400 km it is independent of pitch angle but rises rapidly 
with altitude to the value of 17 eV at 400 km. Above this altitude the 
mean energy is a function of both pitch angle and altitude, rising with 
altitude to about 18 eV at 1000 km for p > 0, while decreasing to about 
14 eV for y < 0. The complicated structure of the mean energy demon- 
strates the importance of calculating the pitch angle distribution of 
the flux, since otherwise this structure cannot be revealed. 

3.4 Integrated Photoelectron Fluxes. 

The integrated photoelectron fluxes *l>"(h), ^^(h), and $^(h) defined 

by equations (4a, b, and c) are shown in Figure 9. The upward-downward 

electron flux anisotropy, noted previously, is clearly seen to persist 

down to about 250 There are three additional conclusions that can 

«»- 

be drawn from the examination of this figure. 

First, the steady-state photoelectron flux attains its maximum 
value at a higher altitude than that where the photoelectron production 
rate is maximum. In section 3.3 it was noted (see Mantas 1973) that 
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the photoelectron production rate for x “ maximum at about 250 km. 

# - 

From Figure 9. it is seen that for x = O'" the total photoelectron flux is 
maximum at about 220 km, v;hile for x = 00® is maximum at about 290 km. 
Therefore, at noon the total flux attains its maximum value about 70 km 
above the altitude where the production rate is maximum. At dawn the 
displacement is 40 I^. This tlisplacement is not due to 

transport. This can be deduced from the observation that the photoelectron 
flux distribution Is isotropic in pitch angle below about 250 km, which . 
implies that there is no net photoelectron current below this altitude, 
and also, from the fact that the displacement is larger in the noon than 
in the dawn model, contrary to what one would expect if the displacement 
were due to transport. The cause of this displacement must, therefore, 
be sought elsewhere. The most likely cause Is an altitude dependent 
photoelectron lifetime. To derive the photoelectron lifetime one would 
have to solve the time dependent transport equation, which at the 
present seems impractical. An indication that the displacement is due 
to this cause ,is given by the altitude dependence of the mean photo- 
electron energies. From Figures 4, 8 and 9 is seen that the photoelectron 
flux maximum occurs at exactly the same altitude as the altitude above 
which the mean photoelectron energy increases rapidly. In Figure 4 
there appear to be two such changes; the altitude of the photoelectron 
flux maximum coincides with that of the upper change. 

Second, the upward photoelectron flux above 500 km, and therefore 
the escape flux, is larger at dawn than at noon. With regard to this 
conclusion it should be noted that the upward flux at 1000 km at noon 
is lower than the maximum upward flux at 220 km by about a factor of 10, while 
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at dawn this factor is about 5. The stronger altitude dependence of 
the upward flux at noon is due to the larger column number density 
(neutral molecules and ions included) at noon than at dawn. 

And* third, the photoelectron flux below 250 km depends strongly 
upon the solar zenith angle, while above this altitude the dependence 
on X Ts weak. Therefore, while comparison between measured and calculated 
fluxes below 250 km is difficult unless the solar zenith angles are the 
same, above 250 km meaningful comparisons between theory and experiment 
can be made even for different solar zenith angles. 

3.5 Photoelectron Escape and Energy Input into the Protonosphere. 

Above 700 km the upward photoel ectron flux $ is seen (Figure 9) to 

change very little with altitude. This implies that very few of the elec 

trons that reach this altitude are backscattered into the ionosphere while 

the majority escape to the protonosphere. To obtain an estimate of the 

backscattered fraction of photoelectrons we consider the rate of change 

with altitude of the ratio ($”/^^) at the base of the protonosphere. In 

the altitude interval 800 to 900 km we find (from Figure 9, x = 90"") that 

-5 -1 

the altitude dependent reflection coefficient is about 5 x 10~ (km“*). 

If we assume that photoelectrons are reflected at the same rate in the 
first 1000 km above the base of the protonosphere, and beyond this 
altitude they are completely absorbed, the reflection coefficient ’r. \ 
at the base of the protonosphere (h = 1000 km) is 0.05. This value of 
R|^ is, at best, a rough estimate; it demonstrates, however, that most 
of the upward flux at high altitudes escapes to the prbtonosphere. In 
the examination of the differential flux spectra (Figures 2 through 2c) 



we noticed that at high aEltltudes the backscattering probability is 

0 

much higher for low than for high energy electrons. This implies that 
the reflection coefficient is energy dependent. The energy dependent 
reflection coefficient can be obtained from the energy spectra by 
arguments similar to those used above. However, our Interest is not 
to calculate the reflection coefficient but simply to show that it is 
small, and proceed with the calculation of the escape flux. 

The net number of photoelectrons per electron volt crossing a unit 
surface perpendicular to the z-axis (the direction of the magnetic field), 
in the direction of increasing z, per second, at the altitude h, is 

^>p(h,E) = 27 t \ u ^(h,E,]j)du (7) 

-1 

where we have inserted the subscript n to distinguish the net flux parallel 

to the magnetic field line defined by' (7) from 0^(h,E) defined by (2). 

Evaluation of (7) at h = 1000 km gives the energy spectrum of the 

escaping flux. A negative flux $^(h,E) at 1000 km implies that a net 

number of electrons enter the ionosphere at the upper boundary. Integration 

of this equation over E at 1000 km gives 2.6 X 10^ cm"^ and 1.5 X 10^ 

sec"^ for the escaping electron fluxes at dawn and noon, respectively. 

Multiplying (7) by E and integrating over E at 1000 km we obtain 4.4 X 10^ 

~2 -1 

eV cm sec for the energy input rate into the protonosphere at dawn, and 
9-2-1 

2.7 X 10 eV cm sec for the energy input rate at noon. These energy 
input rates compare favorably with the downward conduction heat fluxes 
(Evans 1967; Evans and Mantas 1968) obtained from observed electron 

■|*omnor;»tniro nv'aHion+c Af hinh Alti+iirlpc Millcfnnp HtII 
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4. COMPARISON WITH EXPERIMENT 

Before we compare the present calculations with measured photoelectron 
flux energy spectra we like to point out that the theoretical energy spectra 
for X “ shown In the following three figures have been recalculated. 

This was considered necessary in order to evaluate the effect of an error 
in the computer program which was found after most of the calculations had 
been performed. It was found that this error affected only the low energy . 
region of the spectrum > especially in bringing out a more pronounced effect 
on the spectrum due to the vibrational excitation of Since It did not • 
have any substantial effect elsewhere, nor did it affect the conclusions 
drawn from the previous calculations, it was deemed unnecessary to repeat 
all the calculations. 

Before comparing the calculated spectra with experimental results we 
would like to point out the extent to which agreement can be expected. 

First, the energy resolution of the experimental results is not sufficient 
to confirm or disprove the finer structure of the energy sepctra found 
in our calculations. Therefore, find structure comparison must be by- 
passed until experimental data of the required degree of resolution become 
available. Second, the number of experiments is limited, being confined 
to altitudes below 300 km, where the fluxes depend strongly upon the solar 
zenith angle. An exception to this is the data of Galperin ^ (1972) 

which refer to altitudes above 300 km; however. 
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the absolute value of the fluxes measured in these experiments is questionable 
*du6 to calibration problems. These authors have normalized their data to 
the results of Doering et al . (1970), at 300 km, in the energy region 
between 30-40 eV. 

In Figure 10 the calculated photoelectron flux energy spectrum at 

179 km is compared with the experimental data of Doering e^ al^. (1970) 

and Knudsen and Sharp (1972) at 180 km. In the same figure we have also 

included the data by Galperin et al, (1972) at 250 km. The experimental 

_2 -1 -1 -1 

data were given as electrons cm eV ster sec ; here we 
have multiplied these data by 4it to obtain total fluxes. No error is 
introduced by this, since our calculations have shown that the fluxes are 
uniformly distributed over all directions (isotropic) at this altitude. 

The agreement is excellent for energies above 5 eV. Below this energy 
there are only two data points at 4 eV. In the caluclated spectrum the 
slope changes sign at about 5 eV, while in the measured spectrum no such 
change is seen. However, a recent experiment by Hays and Sharp (1573) 
that covers the energy region below 5 eV clearly shows the change in the 
slope between 4 and 5 eV and the presence of a minimum at about 2.5 eV. 

The observed minimum is, however, much shallower than that in the calculated 
spectrum, indicating either insufficient experimental resolution,or an 
overestimated cross section for vibrational excitation of N 2 . It should 
be noted, however, that the overestimate need not be large since only 
small change ir\ the ratio of the inelastic to the elastic cross sections 
affects substantially the calculated spectrum. 

In Figures 10-12 it is interesting 

to note that broad features In the calculated spectrum are also present in 
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,the measured spectra. For example, in the energy region about 15 eV the 
measurements of Doering il- (1970) definitely show structure; this 
is also true In the energy Interval between 40 and 50 eV in the data of 
Knudsen and Sharp (1972). Since in both regions the calculated spectra 
also show broad features, the absence of the sharp peaks between 20 and, 

30 eV in the measured spectra indicates that the experimental resolution - 
is not sufficient to resolve these peaks. 

From Figure 9 it can be seen that the magnitude of the flux at 180 km 
and 250 km is approximately equal (for x = 0°). Therefore,it is not 
inappropriate to include the Cosmos 348 results at 250 km, by Galperin 
(1972), in this comparison. 

From this comparison we can see that the theory is in good agreement 
with experiment at the high energy region also. Between 30 and 60 eV the 
agreement is good. Above this energy the calculated spectrum is lower by 
a' factor of 3, or less, depending on the energy. Above about 70 eV the 
calculated spectrum is changing little with energy, while the measured 
spectrum continues to decrease at a different rate, however, than that 
observed below 70 eV. The change in the slope of the measured spectrum 
rectrum reflects the change seen clearly in the calculated spectrum. 

Taking into account the fact that we have not included energies above 
100 eV, it is expected that the calculated spectrum between 70 and 100 
eV would be lower than the actual distributioD, In this region the 
calculations can be extended straightforward to include energies above 
100 eV. Such an extension would require only calculation of photoelectron 
production rates with energies above 100 eV, and extension of the 
collision cross sections over this energy region. Although an extension 

i 
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that Includes production above 100 eV will affect the energy spectrum 
everywhere, the energy regions that are expected to be mostly affected 
are the high and low energy regions, the first because cascading from 
energies above 100 eV would be nonzero, and the second, mostly through 
secondary low energy electron production by ionization. 

In Figure 11 the calculated spectra are compared with the results 
of Knudsen and Sharp (1972) for several altitudes below 180 km. At 
altitudes above 140 km the agreement is fairly good; however, the change 
in the slope at about 5 eV, found in our calculations, does not appear 
in the data. We have already pointed out reasons for this discrepancy. 
The calculated spectrum at 121 km is lower than the observed. The 
disagreement at this altitude is not too large If we take into considera- 
tion the large solar zenith angle difference, and the fact that the cal- 
culations reprsent low solar activity conditions, while the data corres- 
pond to a higher level of solar activity. 

In Figure 12 we compare calculated spectra at various altitudes vnth 
the measurements of Doering et al_. (1970). Since the large solar zenith 
angle difference precludes detailed comparison, we shall compare only 
the shapes of the calculated and measured spectra. If we exclude the 
change in the sign of the slope of the spectrum at about 5 eV, which as 
we have noted has been observed in a more recent experiment, the agree- 
ment is good. It should be noted here that since the calculations are 
very sensitive to changes in the ratio of the inelastic to the elastic 
cross sections, the agreement shown in this figure can be substantially 
improved by small changes in the cross sections in future calculations. 
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5. CONCLUSIONS 

Theoretical calculations of photoelectron flux energy spectra have 
revealed several important characteristics of the photoelectron distribution 
function, the most important of which are: 

(a) The angular distribution of the photoelectron flux is isotropic 
at all altitudes below about 250 km. Above this altitude an upward- 
downward flux anisotropy is found which increases with altitude. This 
anisotropy gives rise to a larger upward flux component, and therefore 

a net photoelectron current in the upward direction. 

(b) At altitudes above 550 km the mean energy of the upward flux 
component is larger than that of the downward component and changes 

very slowly with altitude. The mean energy of the downward flux component 
on the other hand decreases rapidly with altitude above 550 km. In the 
region between 250 and 550 km the mean energy of the downward flux is 
slightly larger than tftat of the downward component. Below this 
altitude energy, as well as angular, isotropy prevails. 

(c) The total photoelectron flux attains its maximum value at 
altitudes higher than that where the production rate is maximum. Since 
transport has been found insignificant below 250 km this displacement 
has been attributed to an altitude increasing photoelectron lifetime. 

The altitude where the flux is maximum coincides with the altitude above 
which the mean photoelectron energy begins to rise rapidly with altitude. 
Therefore, photoelectrons at altitudes higher than that where the flux 

is maximum are on the average more energetic than those below this altitude. 
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(d) The total number of photoelectrons escaping to the protonosphere 
at dawn is larger than the number escaping at noon by about a factor of 
two. The same is also true with regard to the energy supplied to the 
protonosphere by these electrons. 
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FIGURES 

• 

Figure 1. Energy spectrum of the primary photoelectrons at selected 
altitudes at dawn [y = 90°] {adapted from calculations by 
Stolarski 1972, private communication). 

Figures 2a - 2c. Photoelectron pitch angle and energy spectra at 

selected altitudes at dawn. For negative values of p ly = cosine 
(pitch angle)] the electrons are moving In the downward direction, 
while for positive values they are moving upward. Note the pitch 
„ angle isotropy at low altitudes (a,b) and upward -downward uniso- 

tropy at high altitudes. Also, note the broadening and shifting 
towards lower energies of the spectrum features around 25 eV 
with decreasing |p| at high altitudes and, also, the difference 
in the sign of slopes of the upward and downward spectra at the 
low energy limit. 

Figure 3. Photoelectron flux energy spectrum at dawn. Note the gradual 
smoothening of the spectrum with increasing altitude and, also, 
the broadening and shifting towards lower energies of the 
spectrum features around 25 eV. 

Figure 4. Mean photoelectron flux energies: is the mean energy of ■ 

the primary photoelectrons (source), top scale is the mean 
energy of the downward photoelectron flux, and V‘ that of the 
upward, bottom scale. The curve to the right represents the 
ratio E*^/E-. . The mean photoelectron flux energy follows 
closely the F** curve. Note that while the mean energy of the 
primary photoelectrohs increases with decreasing altitude, the 
mean energy of the steady-state photoelectron flux decreases. 
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Figure 5. 


Figure 6. 


Figure 7. 


Figure 8. 


Figure 9. 


Figure 10, 


Figure 11. 


Figure 12, 


Pitch angle distribution of the total (integrated over all 
energies) photoelectron flux. Note the large upward -downward 
anisotropy jLbove 500 km. 

Photoelectron flux energy spectrum at noon. Note the rapid 
transition in the smoothness of the spectrum at about 300 km, 
and the broadening and shifting of the spectrum features with 
increasing altitude {compare with Figure 3), 

Energy spectra of the upward flux at selected altitudes at 
noon. The shift and broadening of the features at 22,5 eV, 
and the energy dependent attenuation of the flux with 
increasing altitude, are clearly seen in this figure. 

Mean energy of the steady-state photoelectron flux at noon 
as function of altitude and pitch angle cosine. 

Total photoelectron fluxes at dawn and noon. Note the upv/ard- 
downward flux unisotropy above about 250 km. 

Comparison between calculated and measured photoelectron flux 
energy spectra at 180 km. 

Comparison between calculated and measured photoelectron flux 
energy spectra at selected altitudes between 120 and 180 km. 
Comparison between calculated and measured photoelectron 
flux energy spectra at selected altitudes between 120 and 303 km. 
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